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We extend a method for modeling synchrotron and synchrotron self-Compton radiations in blazar 
jets to include external Compton processes. The basic model assumption is that the blazar radio 
through soft X-ray flux is nonthermal synchrotron radiation emitted by isotropically-distributed elec- 
trons in the randomly directed magnetic field of outflowing relativistic blazar jet plasma. Thus the 
electron distribution is given by the synchrotron spectrum, depending only on the Doppler factor Sr> 
and mean magnetic field B, given that the comoving emission region size scale R' b < c5T>t v /(l + z), 
where t v is variability time and z is source redshift. Generalizing the approach of Georganopoulos, 
Kirk, and Mastichiadis (2001) to arbitrary anisotropic target radiation fields, we use the electron 
spectrum implied by the synchrotron component to derive accurate Compton-scattered 7-ray spectra 
throughout the Thomson and Klein-Nishina regimes for external Compton scattering processes. We 
derive and calculate accurate 7-ray spectra produced by relativistic electrons that Compton-scatter (i) 
a point source of radiation located radially behind the jet, (ii) photons from a thermal Shakura-Sunyaev 
accretion disk and (iii) target photons from the central source scattered by a spherically-symmetric 
shell of broad line region (BLR) gas. Calculations of broadband spectral energy distributions from 
the radio through 7-ray regimes are presented, which include self-consistent 77 absorption on the 
same radiation fields that provide target photons for Compton scattering. Application of this base- 
line flat spectrum radio/7-ray quasar model is considered in view of data from 7-ray telescopes and 
, contemporaneous multi-wavelength campaigns. 

Subject headings: radiation processes: nonthermal — galaxies: active — supermassive black holes 
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1. INTRODUCTION 

A class of radio- loud blazar active galactic nuclei (AGNs) that emit luminous fluxes of > 100 MeV - GeV 7 



s . rays was discovered with the Energetic Gamma Ray Experim ent Telescope (EGRET) on the Compton Observatory 
(jHartman et al.lll992t iFichtel et alJll994t lHartman et~aOl999| ). Thi s result clarified the n ature of 3C 273, which was 
f - ^ ' first identified as a 7-ray emitting AGN in COS-B satellite data (Hermsen eTa!][l977h . The 7 rays from blazars 
are certainly nonthermal in origin and associated with the radio jets formed by the supermassive black holes that 
power these sources. The largest subclass of EGRET AGNs are moderate redshift (z w 1) flat-spectrum radio quasars 
(FSRQs) with blazar properties, including apparent superluminal motion, rapidly variable optical emission, high 
polarization, and intense broadened optical emission lines. Another subclass of the EGRET AGNs consists of 7-ray 
' emitting BL Lacertae (BL Lac) objects, which are generally at lower redshifts (z w 0.1 - 0.3) and, by definition, have 
weak or absent optical emission lines in their spectra. The X-ray selected BL Lac (XBL) su bset were discovered to be 
a class of TeV 7-ray sources by the Whipple Observatory (|Punch et al.lll99 2: Weekes 20Ql). 

The common superlum inal nature of the first identified 7-ray blazars, namely 3C 273, 3C 279, and PKS 0528+134, 
led lDermer et al.l {l992) to propose a Compton-scattering origin for the 7 rays. In this model, jet electrons Compton- 
scatter accretion-disk photons that intercept the jet plasma. The nonthermal j et electrons can also scat ter internal 
synchrotron photons to produce a synchrotron self-Compton (SSC) component (jBloom fc Marscherlll996l ). Given the 
broadened emission lines in the spectra of FSRQs, accretion-disk radiation scattered by surroundi ng gas of the broa d 
line region (BLR) will provide a further source of t arget photons to be scattered to 7-ray energie s (jSikora et al.l ll994). 
as will radiation from a surrounding dusty torus (|Kataoka et a l. 1999; Bla zeiowski et al.ll2000f). The accretion-disk 



and scattered radiation wil l attenuate jet 7-rays through 77 pair-production attenuation (Becker & Kafatos 1995; 
iBlandford fc Levinsonlfl99l . 

Expressions for the 7-ray spectral energy distributions (SEDs) of blazars produced by Compton scattering processes 
have been derived and calculated for many specific models of the black- hole/blazar jet environment. In the case of 
external accretion-disk photons as the target photon source, where the accretion disk is described by an optically 
thick, geometrically thin thermal IShakura and Sunvaevl (I1973T) accre t ion d isk, Compton-scattered 7-ray spectra were 
calculated in the Thomson regime by iDermer fc Schlickciscr (1993], l2002f ) . Calculations of the Thomson-scattere d 
spectra for a quasi-isotropic target radiation field formed by BLR gas or hot dust were made by ISikora et al.l (|1994f> . 
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iDermer et al.l ((1997T ) and iBIazeiowski et all ([20001 ) . D etailed numerical calculation s including both accretion disk 
and scattered radia t ion fi elds, have been made by, e.g., iKusunose fe Takaharal (|2005l ): lB5ttcher fc Blooml (|2000l ) and 
iBottcher fc Reimerl (|2004 . 

Compton scattering in the Klein-Nishina regime is not so simple to treat compared to analyses restricted to the 
Thomson regime, but is unavoidable for blazar analysis in the era of the Fermi Gamma-Ray Space Telescope (formerly 
known as GLAST) and the ground-ba sed 7 Cherenkov telescopes. For surrounding isotropic radiation fields in the 
stationary frame of the blazar AGN, IGeorganopoulos et alj (|2001l) suggested to transform the comoving electron 
distribution to th e stationary frame and then scatter th e target photons to 7-ray energies, using the formula first 
derived by Jones (|Jones 1968; Blumenthal fc G ould 1970f). This approach is generalized in this paper to surrounding 
anisotropic radiation fields. 

The usual spectral modeling approach proceeds by injecting power-law electrons and evolving these particles 
while they produce the output syn chrotron and Compto n-scattered radiation (e.g., IDermer fc Schlickeiserl 119931 : 
lB5ttcher et al.l fl997f ) . For example, iModerski et al.l (|2005f ) calculate electron energy evolution and spectral forma- 
tion throughout the Thomson and Klein-Nishina regimes for different ratios of synchrotron and isotropic radiation 
field energy densities. They show that reduced Compton losses in the Klein-Nishina regime compared to synchrotron 
losses can lead to spect ral hardening of the synchrotron component in the optical/X-ray regime (noted earlier by 
IDermer fc Atovanll2002D . A difficulty in this approach is that the electron energy-loss rate depends on the photon 
spectrum of the comoving radiation field, not just the total radiation energy density, and this field evolves with time. 
The modeler is faced with the prospect of simultaneously fitting the synchrotron and Compton components. The 
acceleration scenario may well be over-simplified, and non-power-law particle injection distributions could be more 
realistic than power-law injection spectra, e.g., due to nonlinear effects in Fermi acceleration. Moreover, a separation 
between the acceleration and radiation zones may not be justified. 

Here we extend a method of blazar analysis recently proposed for TeV blazars (jFinke et al.l [2008) that avoids these 
difficulties. For a standard 7-ray blazar model, where isotropically distributed electrons spiral in a randomly oriented 
magnetic field with mean magnetic field strength B in the fluid frame, the measured synchrotron flux directly reveals the 
electron spectrum responsible for the synchrotron radiation. The only uncertainties are the mean magnetic field B, the 
comoving size scale R' b of the emitting region, and the Doppler factor 5^ = — j3 cos 9)]^ 1 (T = \j\J\ — ft 2 is the bulk 
Lorentz factor of the outflow). With this electron spectrum, we then Compton-scatt er target photons of the surrou nding 
radiation fields using the head-on approximation to the total Compton cro ss section (jDermer fc Schlickciscr 1993), valid 
when the electron Lorentz factor 7> 1. This generalizes the approach of IGeorganopoulos et alJ (|2001l ) to surrounding 
anisotropic radiation fields. The temporally-evolving electron spectrum in blazars can be derived in this approach 
from simultaneous multiwavelength blazar data. Values of B, 5b, and jet power can then be deduced. The related 
treatment for XBLs applied to PKS 2155-304, including more details about the derivation of the electron spectrum 
from the synchrotron component, th e derivation and ca lculations of the SSC component and internal 7-ray opacity by 
the synchrotron photons, is given by iFinke et al.l (|2008f ) . 

Analysis of blazar SEDs using this approach is presented in Section 2, where formulas to calculate Compton-scattered 
internal and external radiation and a ^-function approximation for 77 opacity from the internal radiation field are 
given. Derivations of the Compton-scattered spectrum for specific examples of external radiation fields consisting of 
a monochromatic point source of radiation radially behind the jet, a Shakura-Sunyaev disk model, and a model BLR 
radiation field are derived in Section 3. Discussion of the results is found in Section 4. 

2. ONE-ZONE SYNCHROTRON/SYNCHROTRON SELF-COMPTON MODEL WITH 77 OPACITY 

We consider a one-zone model for blazar flares. Multiple zones could still be allowed, but the product of the duty 
cycle and number of zones would have to be small enough that interference of emissions from the different regions would 
still permit rapid variability. In this case, the emission would still predominantly arise from a single zone. Distinct 
zones could also emit the bulk of their radiation in different wavebands. In this case, the cospatiality assumption often 
made in blazar modeling would not apply. In this regard, correlated variability data is essential to test the underlying 
assumptions made when a one-zone model is employed. Slowly varying radio/IR synchrotron and hard X-ray and 
low-energy 7-ray Compton emissions could involve extended emission regions. 

A radiative event from the source emission region that varies on a comoving timescale t' v > R^/c is related to the 
observed variability timescale through the relation t' v — 6jyt v /(l + z), where z is redshift; thus the comoving blob 
radius is R' h < c5r>t v /(l + z). The inequality allows us to neglect light-travel time effects from different parts of the 
emitting volume and avoid integrations over source volume. Within this zone, the nonthermal electrons with isotropic 
pitch angle distribution are described by the total comoving electron number spectrum N^j') in terms of comoving 
electron Lorentz factor 7'. The magnetic field is assumed to be randomly oriented in the comoving fluid frame. The 
relativistic electrons that gyrate in this field radiate nonthermal synchrotron radiation, observed as the low-energy 
component in blazar SEDs. 

2.1. Synchrotron and Self-Compton Components 
The vF v synchrotron radiation spectrum can be approximated by the expression 

f syn ^ _A_ caTUBi m N , M ^ (1) 



7-ray FSRQ Blazar Analysis 



3 



where 



di = (1l(z) is the luminosity distance, c is the speed of light, ot is the Thomson cross-section, z is the source redshift, 
and the comoving magnetic-field energy density of the randomly-oriented comoving field with comoving mean intensity 
B is Ub = B 2 /8ir . We use e and e' to refer to the dimensionless photon energy in the observer and comoving frame, 
respectively. Here and throughout this paper, unprimed quantities refer to the observer's frame, and primed quantities 
refer to the frame comoving with the AGN's jet, with the exception being B, the comoving magnetic field. Inverting 
this expression gives the comoving electron distribution 



6nd 2 L f, svn 
c<jtUb$vT. 

where 



n'M = HW) = —r, AZa > ( 3 ) 



e S yn 1 + z , W 

cb = B/B cr is the ratio of B and the critical magnetic field B cr = m 2 c 3 /eh = 4.41 x 10 13 G (jDermer fe Schlickeiserl 
120021 ) . and V£ = is the comoving volume of the blob. Note that U B = e 2 B U Bcr = e 2 B B 2 r /8ir. Eq. © gives a 

good representation to the source electron distribution when the vF v spectral index a < 4/3 (i.e., for spectra softer 
tha n a = 4/3, adopti ng the convention f e cx e a ) and away from the high-energy exponential cutoff of the spectrum 
(see lFinke et al.ll2008l for comparison). 
The SSC vF v flux is given by 

/f 5C = | <L' ssc (e' s ,Q's) ■ (5) 
The formula of lJonei (fl968l ) (see also iBlumenthal fc Gouldl (fl970h ) gives the SSC vF v flux, 

where 



2 q ']nq> + (1 + 2</)(l - q') + \ (^V) (1 ~ 90 



47 



e' Ay' 

and r' e = 4e'7' . (8) 



F e (l-ei/7') 

The synchrotron photons provide a target radiation field with spectral energy density 

u'(e') = eWH^O = , (9) 

using eq. ([2]). The scattered photon energy in the comoving frame is related to the observed photon energy hv = m e c 2 e 
by the relation 

4 = il±il£,f. do) 

From the limits on the integration over 7 ' implied by the limits on q' we find 

(11) 

A \ V fc " s ) 

and 

e' e ' 

Tmaz = , _ S / g(e - e s ) + 72^(4 - e ') (12) 
e t s 

(see iFinke et al.l 120081 for a detailed derivation of synchrotron/SSC models and application to blazars). Here the 
maximum lepton Lorentz factor injected into the radiating fluid is 7 2 ', and the Heaviside function H (x; a, b) is defined 
such that H (x; a, b) — 1 when a < x < b, and iJ (x; a, 6) = otherwise; the Heaviside function with one entry is defined 
such that H(x) = 1 when x > 0, and H(x) — otherwise. The vF v SSC spectrum is therefore given by 

fSSC (z\ 3 die' 2 r,,rr , w.ra/r n ^ 
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where e = 8-ne'/{l + z) and e = Sd€bJ i2 /(1 + z). The maximum vF v SSC flux at photon energy e^ c can be 
approximated in the Thomson limit by the expression 

rSSC „ 247Td|(l+z) 2 / syn \ 2 

% c - {8lBt v y<? V<T) ' (14) 

where the peak frequencies are related by 



c 



& ~ V J TT^ (15) 

(|Tavecchio et al.lll998t iFinke et al.ll2008f ). Here f^FJt is the i/i*!, peak of the synchrotron component, which reaches its 

, syn 

maximum at e = e pk . 

Second-order SSC takes place when the SSC photons are again Compton scattered by electrons in the same blob, an d 
may account for superquadratic variability of the 7-ray flux with respect to the synchrotron flux (jPerlman et alj |2008) . 
In principle, these photons can again be Compton scattered to arbitrarily higher orders, though higher-order scatterings 
are negligible due to Klein-Nishina effects. Calculating second-order SSC can be done accurately by replacing 

ff yn with fi sc in eq. (QD, so that 



f ssca = ( n 3 die? r- g?c Fc{q >x e )fr (16) 



2.2. 77 Opacity 

Gamma-ray photons are subject to 77 attenuation by synchrotron photons produced in the radiating plasma, by 
ambient photons in the environment of the black hole (starred frame), and by photons of the intergalactic radiation 
field. The 77 pair-production cross section 



(17) 



(jJauch k Rohrlich|[l976l : lNikishovlll96ll ; iGould fc Sch7idenfl967t 1 rown. Mikaelian. &: Gouldlfl973t ). where 7 cm is the 
center-of-momentum frame Lorentz factor of the produced electron and positron, (3 cm = (1 — 7„ 2 ) 1 ^ 2 = Vl — s _1 , 



, _ e*ei(l + z) 



s = 7cm = ^ ^ ~~ C ° S ^^ ' ^ 

and r e = e 2 /m e c 2 = 2.8179 x 10 -13 cm is the classical electron radius. The interaction angle i/j, given by the relation 



COS?/' = fX^ s + y/l - fll\Jl - /1 2 COS(0» - <t>s) , (19) 

is the angle between the directions of the photon detected with energy e\ and the target photon with energy e* . 
The absorption probability per unit pathlength is 

— 22 — — = j> dQ* (1 - cosf/ 1 ) J de* n ph (e*,£l r ) cr T7 (s) . (20) 

For absorption by synchrotron photons within the radiating volume, e* — > e' and ex — > e[ = (l + z^/Jrj, and the target 
synchrotron radiation held is given by eq. In this case, the optically-thin 7-ray emission spectrum is modified by 
the factor 3w(t 77 )/t 77 for a spherical geometry, where 



1 exp(-r 77 ) 1 - exp(-r 77 ) 

U ( T 77) = ^ " ~2 ' ( 21 ) 

z '77 '77 

Here t 77 is the total 77 optical depth integrated over pathlength. For absorption by ambient photons in the vicinity of 
the AGN, is the photon energy in the AGN rest frame. For cosmic 77 absorption, the target photons are given by 
the spectrum of the intergalactic background light, which evolves with redshift. In the latter two cases, the intrinsic 
spectrum is modified by the factor exp(— t 77 ). 

3. COMPTON-SCATTERED EXTERNAL RADIATION FIELDS 

In the one-zone model, the vF v spectrum of Compton-scattered external radiation fields is given by the Compton 
spectral luminosity e s Lc(e s ,^s) according to the relation 

f? = £sL °| S ' ns) , (22) 

where e s = (1 + z)e, from eq. (fit)]) , and fl s = il. The latter equality means that the photon direction is not deflected 
in transit to the observer. The Compton spectral luminosity is given by 

/> />oo /• />oo 

e s L c (e a ,Q s ) = m e c 3 e 2 <b dtl* / de* n ph (e*, Q*) * dQ e / d-y N e (j,Q e )x 
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(i_ cos ^)^°(f)$(n,-n e ) ) (23) 

de s 

having already introduced the approximation that the scattered photon travels in the same direction as the relativistic 
scattering electron, i.e., Q, s = il e . Because of this approximation, the cosine of the angle ip is given by eq. Q19p . The 
invariant collision energy 

e = 7e*(l — yl — l/7 2 cosf/') — 7£*(1 — cosip) (24) 

because 7 > 1. The relation n p /j(e»,il*) = u(e*, f2*)/(m e c 2 e*) gives the specific spectral number density of target 
photons with energy e*, the starred quantities referring to the frame stationary with respect to the black hole. 
The Compton cross section in the head-on approximation is given by 

da c ^ irrl ( e 2 7 e \ 

(jDermer fc Schlickeiserlll993l : iDermer fc Bottcherll2006l ). where 

S = y + y-i-^ + ^a, (26) 
7ey 7ey 

y = 1 - - , (27) 

7 

e is given by eq. (I24[) . The Compton spectral luminosity in the head-on approximation becomes 

dQ* / &. * l * 2 ' * y / d 7 7-^(7, fi s )S. (28) 
The lower limit on the electron Lorentz factor 7; ou , and the upper limit e^M implied by the kinematic limits on y are 

(29) 



"flow — 2 

and 



1 + 4 1 



e*e s (l — costp) 



e*,W = — — — r ■ (30) 
1 — cos %p 

Eq. (P5|) is the starting point to calculate accurate Compton-scattered spectra involving relativistic electrons and 
external photon fields with arbitrary anisotropics. In contrast to the comoving electron spectrum used in the SSC 
calculation, the calculation of Compton-sc attered radiation uses the elec tron spectrum N e (j, fi e ) and the target photon 
spectrum defined in the stationary frame (jGeorganopoulos et al.ll200lh . The invariant phase volume dV = dVd 3 p for 
relativistic particles is given by 

dN _ dN 1 1 dN 

~av ~ dv¥f ~ (m e c) 3 7^ djdndv ~ mv ' ( ' 

implying that 

N ^ Q ) = Ktftf) = SlKb'M) , (32) 

noting that dV/dV — dt'/dt = <5d, and 7 = 5d7' when 7', 7 S> 1, required for the head-on approximation. For an 
isotropic comoving distribution of electrons, iV e (7, f2 s ) = S^N^j')/^. Hence 

e s L c (e s ,n s ) = °4- &l #* t f "'" ** ^%^> j°° d 1 7 ~ 2 ^(7/<5d) 3 , (33) 



* Jo J-l Jo 

or 

fEC _ e s L G (e s ,n s ) _ cirri 2,3 A , f\ F M , f° , Kb/fa) „ 



>. = ^f^ = ^H fc^jf^^B. (34) 

In terms of the measured synchrotron i/J^ spectrum, eq. © , the source Compton spectrum for external Compton 
(EC) scattering in a standard one-zone model for blazars is, in general, given by the four-fold integral 



with 

2 

g = pi 5 u ' ( 3fi ) 

(1 + Z)0d 

using eq. ([3]). The number of integrations can obviously be reduced by choosing symmetrical target photon geometries. 



6 



Dermer, Finke, Krug, and Bottcher 



3.1. Point Source Radially Behind Jet 

First we consider the flux when nonthermal electrons in a relativistic jet Compton-scatter photons from a point 
source of radiation, isotropically emitting and located radially behind the outflowing plasma jet. For a monochromatic 
point source with luminosity L and energy e , the spectral luminosity can be expressed as 

£*(e*) = L Q 5(e* - e ) . (37) 

The spectral energy distribution of the target photon source at distance r from the point source is therefore given by 

«(e., n.) = tb 1 *- 1 ) ^ } . (3g) 

47rr z c Z7r 

Substituting eq. (f3"5| into eq. (j53"|) and solving gives 

a5(( ,, n , )= !|Mr, T fflWi. (39) 



eC, P t _ 



M (40) 



16nr 2 d 2 L eg 7 2 



Using eq. (fl2"2")) . eq. (|3T))) becomes 
or with eq. 

where e is defined in eq. (f3"r?)) . S is defined by eq. (|2rJ)) with e replaced by e = 760 (1 — Ms), an d 



/^ = lz^f% / (41) 



^low ' 2 




e e s (l - Ms) 



(42) 



Eqs. (j40|) and (jITj) give the Compton-scatt ered spectrum from a point source of radiation located radially behind the 
jet, generalizing the Thomson-regime result (jDermer et al.lfl992t ) to include scattering in the Klein-Nishina regime. A 
scattered disk component should be found in all blazar models, with its importance strongly dependent on distance r 
of the jet from the accretion disk. The point-source approximation gives the least upscattered flux in the Thomson 
limit, and an extended disk having the same power as a point source will give a more intense flux. At sufficiently large 
jet heights r > T 4 i? g , defining the far field, where R g = GM/c 2 = 1.5 x 10 13 M 8 cm is the gravitational radius, the 
Shakura-Sunyaev disk can be described as a point source radially beh ind the the jet. Photons from large disk radii 
are important in the near field r < T 4 R g CE ermer ASchlickciscr'20Ql). 

3.1.1. Reduction to the Thomson Regime 

We now derive the Thomson limit for the vF„ spectrum, eq. (j40|) . Because we consider relativistic electrons 7 ^> 1, we 
are restricted to the condition ji ow ^> 1, which occurs according to eq. (|4*2"1) when either e s ^> 1 or e s /eo(l— fJ, s ) ^S> 1. The 
Thomson condition can be expressed as e s <C 7, which is guaranteed when e s <C ji ow , in which case 2eoe s (l — fj, s ) <C 1. 
Another statement of the Thomson condition is that 760(1 — /x s ) <C 1 which, with e s -C 7, again implies that 
e e s (l-M s ) < !• Thus 



How -> x 7r~n t • ( 43 ) 

2e (l - Us) 



For the scattering kernel, eq. (|26p , f s « 7 and y — > 1 in the Thomson regime, so 



s7 e / \7 e 

Away from the endpoints of the spectrum, e s <C 7c and St — > 2. Hence 



S->St = 2-2(S) + (5] . (44) 



/f ,T 

defining 



^v*r*r^. (45) 



4 (imi L y V«o/ i&i T 2 



For the comoving electron distribution, eq. ([3]), in the power-law form 

N' e { 1 ')=K'i^H{ 1 '- 1 i 1 i) 1 (47) 
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eq. (|45|) becomes 



4(p + l) (47rrd L ) 2 V^o 



2 



'D 



max(7,7 1 ) - 7z ^ 



a T L K' ( e s 



/ c \ (3-p)/2 

(|J C P [2(1-M)] (P+1)/2 , (48) 



4(p+l) (47rrd L ) 2 V e o 
where the final expression applies in the regime 7/ < 7 < 72 • This can be written as 



which can be compared with the Thomson-regime expression 
(|Dermer et alj|1992t iDermer fc Schlickeisenfl993l I2002D . where 



7^(7') (50) 



~ / 



1 / e(l + z) 



Sd V £ o(l - fJ-s) 



= V2 7 . (51) 



3.1.2. Accurate Thomson Regime Spectrum 

Eq. (|49|) for the Thomson-scattered spectrum was derived assuming St = 2, away from the endpoints of the spectrum. 
Using the full expression for St and a power-law electron distribution gives an accurate expression for the spectrum 
of a localized jet of isotropically entrained electrons Thomson scattering a point source, monochromatic radiation field 
that enters the jet from behind. The result is 

where 

S = e (1 S a ) = 2&lf ■ (53) 
£o(4 - fi s ) 

For the power-law electron distribution N' e (j'), eq. (|4"T)) . the accurate analytic Thomson- regime vF v flux from an 
isotropic monochromatic point source radiation field located behind the jet is 



f pt.T _ 3 (ttLqK' .c,.3 ,3+p f 1 r -(1+p) _ - 



- 7 2 - (3+P) ] + ^-r [7r (5+P) - 7 2 " (5+P) ] }, (54) 



(3+p) L,i " J 2(5 +p) 

where 71 = max(<5D7i, <5d7) and 72 = <5d72'- 
In the asymptotic limit 7l ' <§; 7 <§; 72 ', eq. (j5~4")) approaches 



/f ' T "> ^ (4^) ^(1 " ^ f^'(7) , (55) 

where 



.A(p) ^3 - + . (56) 

\l+p 6+p 5+pJ 

The values of A(p) — 1.0, 0.657, and 0.5 for p = 1,2, and 3, respectively. For the Thomson approximation away from 
the endpoints of the spectrum, given by eq. l[4*9")) . the corresponding coefficient is 3/{p+ 1). 

3.1.3. Solution with Compton Cross Section 

From eqs. §9$) and J22J), 



(47rrd L ) 2 Veo/ 4„„ 7 2 



y + J /- i -^ + (^) 2 

7ey 7 ey 



(57) 



Introducing u — e s /7 and i> = e s en(l — fj, s ) and changing variables to u gives, for the power-law electron distribution, 
eq. 63), 

/e " (4^)2 ^2 J ^ C « IC ' ^ 



where 



c 



du 
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u p 2uP+ 2 



u p - u p+1 



,P+4 



and u\ — £ s 5d/72j u 2 = £s min 



( 1 5d 

V7w ' t7 



,p+2 



1 — it v(X — u) v 2 (l — u)'- 
The series solution of eq. (|59[) is given by 

i+ P +1 2it i+f,+3 1 (i + l)u l+ P+ 5 



P 



1 



E 

i=0 



i + p + 1 v(i + p + 3) 



i + p + 5 



(59) 



(60) 



Eq. (|59|) can be solved analytically for integral p. 



3.2. Shakura-Sunyaev Accretion Disk Field 

For the emission spectrum of an accretion disk surroun d ing th e supermassive black hole, we consider the cool, 
optically-thick blackbody solution of IShakura and Sunvaevl (|1973[ ). The disk emission is approximated by a surface 
radiating at the blackbody temperature associated with the local energy dissipation rate per unit surface area, which 
is derived from considerations of viscous dissipation of the gravitational potential energy of the accreting material. 
The accretion luminosity is defined in terms of the Eddington ratio 



^Edd 



r\mc 



Lb. 



dd 



(61) 



where rj ~ 0.1 is the efficiency to transform accreted matter to escaping radiant energy. The Eddington luminosity 
^Edd = 1-26 x 10 46 Mg ergs s _1 , where the mass of the central supermassive black hole is M — 1Q s MsMq and the 
black hole is accreting mass at the rate m (gm s _1 ). 

For steady flows where the energy is derived from the viscous dissipation of the gravitational potential energy of the 
accreting matter, the radiant surface-energy flux 



d£ 



3GMr 



<p(R) 



(62) 



dAdt 8ttR 3 
(|Shakura and Sunvaevlll973[) . where 

v{R) = [l-p i {R i /R) 1 ' 2 ] , (63) 

Pi = 1, and Ri — 6GM/c 2 for the Schwarzschild metric. Integrating equation (|62[) over a two-sided disk gives r\ = 1/12. 
Assuming that the disk is an optically-thick blackbody, the effective temperature of the disk can be determined by 
equating equation (|62|) with the surface energy flux a^T^[K). A monochromatic approximation for the mean photon 
energy m e c 2 e(R) = fc B T(i?) at radius R of the accretion disk with mean temperature T(R) is given by 



m e c 2 (e(i?)> ~ 2.70k B T(R) 



2.70fc B [ 3GM T ( ^] 1/4 =137(^) 1/4 ^/4 



eV 



(e(R)) ^ 2.7 x 10" 4 £R 



-3/4 



where 



(t 



Edd 



1/4 



VMM, 

Here 9 = arccos/z* is the angle between the directions of the jet and the photon that intercepts the jet, 

R 

Ra 



R 



- 1 



(64) 
(65) 

(66) 
(67) 



where f — r/R g . Lengths marked with a tilde are normalized to R g . The final expression in eq. (f64|) is valid when 

R 3> R m in , though we assume it is reasonably accurate to R = R m in ■ 
The intensity of the Shakura-Sunyaev accretion disk model along the jet axis is given by 

2.7fc B 



(|Shapiro fc Teukolskvl[l983l) . Thus 



rss/^ r,s ^ 3GMm ,_. rr 

3 ^Edd-^Edd 



m e c^ 



T(R)} 



lf s (Sl;R) 



<p(R)S[e - (e(i?)}] 



16tt 2 V RR 2 

(|Dermer fc Sc hlickciscr 2002). Substituting eq. (f69|) into eq. ([33]) . using the relation I £ (£l) — cu(e, Q), gives 



(68) 



(69) 



ss 



er T et 



2 9 tt 3 d\R 2 g 



fedd^Edd 

yyf 3 



ft 



d(f>* 



f{R) 



{^ 2 - lf/*{e{R)Y 
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r d 7 T - 2 ^( 7 /fc) [y + y- 1 - ^ + (4^) 2 ] • 



■yey jcy 

The integration over angle in eq. (|70[) is limited by the inner radius of the accretion disk, so that /i* 

N21-V2 



(6/f) 2 



and 



for a Schwarzschild black hole, and 

I = 1(7, e s ,ip) =7(e(-R))(l-cos^>) , 




(e(R))e s (l-costp) / 

The other limit on the angular integration arises because of the restriction given by eq. (|30[) . so that 

(e(R)) < 2 '' 



(70) 
= [1 + 

(71) 
(72) 

(73) 



1 — cos ijj 

which restricts the integral to a maximum value of R and therefore /i* > l^*. m in- In the calculation of cos?/;, eq. (|19p . 
we take <j) s = without loss of generality because of the assumed azimuthal symmetry of the accretion-disk emission. 

The result for the accretion-disk radiation field scattered by isotropic, relativistic jet electrons is a 3-fold integral — 
reduced from a 4-fold integral by approximating the disk blackbody spectrum by its mean thermal energy at different 
radii. When expressed in terms of the measured synchrotron vF v spectrum using eq. ((3]), the result for the accretion- 
disk radiation field is 



ss 



2 8 tt 2 cR 2 g U B 



fedd-^Edd j-2 
^5 "B 



2 71 



(if i 



<p(R) 



(,i^ - l^eOR))* 



-1 i ( €s \2 

y + y — — + (— ) 

jey jey _ 



(74) 



recalling the definitions of 7J and e from eqs. ([2]) and (f36|) . respectively. This can be reduced to a 2- fold integral by 
approximating a typical scattering as occurring at 0» — n/2, so that e = 7e(l — /i*/i s ). Because it is feasible to perform 
the 3-fold integral numerically, we show results of the more accurate calculations. 

3.2.1. Regimes in Compton- Scattered Accretion-Disk Spectra 

The limiting behaviors of the Compton scattered spectra can be understood based on simple ^-approximations in the 
Thomson regime. External Compton scattering of the Shakura-Sunyaev disk with radiant luminosity Ld = ^Edd-^Edd 
in the near field, i.e., f <C T 4 , can be approximated as 

s2 T 

ECNF ~ "p-^ri 



f, 



ft 



syn 



^syn 



where 



In the far field (r > T 4 



where 



(|Dermer fc Schlickeiseij[2002T) . 



syn(f) 



2ee B 



<5 D e(V3r) 



pECFF 



3 x 10 



-10 



eB(G) (Mtfft 



,~3 \ V4 



Edd 



3 L d 
8 r 2 cB 2 Sl 



3.2 x 10 14 S D ee B M t 



i' Edd 



1/4 



(75) 



(76) 



(77) 



(78) 



3.2.2. 77 Opacity from Accretion Disk Photons 

Photons from the accretion disk will interact with high energy 7-rays to produce electron-positron pairs, modifying 
the very-high energy (VHE; multi-GeV - TeV) 7-ray spectrum by a factor of e~ T7 "'. The absorption opacity r 77 can 
be calculated by inserting the photon density, n p h for an accretion disk into eq. (|20p and integrating x from r to 00. 
For a Shakura-Sunyaev accretion disk, the photon density is given by 



.ss 



(6* , £7* 



e*m e c° 



(79) 
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TABLE 1 

Parameters of Baseline Model 



Rcdshitt 


z 


l 


Bulk Lorentz Factor 


r 


25 


Doppler Factor 


So 


25 


Magnetic Field 


13 


1 G 


Variability Timescale 


t v 


10 4 s 


Black Hole Mass (1O 8 M ) 


M a 


1 


Eddington Ratio 


I em 


1 


Jet Height (units of R g ) 


f 


10 3 


Low-Energy Electron Spectral Index 


Pi 


2 


High-Energy Electron Spectral Index 


VI 


4 


Minimum Electron Lorentz Factor 


7' - 


10 2 


Maximum Electron Lorentz Factor 


n fmax 


10 7 


Accretion Efficiency 


V 


1/12 


Ratio of B to Equipartition Field 


B/B eq 


0.3 


Jet Power in Magnetic Field 




7 X 10 43 erg s" 1 


Jet Power in Particles 


P. 

1 j,par 


10 46 erg s" 1 



where lf s (Q;R) is given by eq. (|69|) . For an azimuthally symmetric accretion disk, the optical depth to 77 pair 
production attenuation for a photon with observed energy ei traveling outward along the jet axis starting at height r 
is given by 



;3/4 „l/4 

r£f(e^) = 3x 10 6 Edd ^ , , , . 
77 V ' ?? 3 / 4 Jf % 2 J 6 R 5/4 (l + R 2 /x 2 ) 3 / 2 



dx f°° dR [^(R)] 1 ^ (s - 1) 



y(5) 



(l-/i*), (80) 



where s = (e(i2))ei(l + z)(l — /tt*)/2 and /i* = l/^J 1 + R 2 /x 2 . For the Shakura-Sunyaev disk extending to the inner- 
most stable orbit of a Schwarzschild black hole, one sees from eq. (|80|) and the definition of (e(R)), eq. ([65]) . and £, eq. 
p)j) , that 

-^1 . — oc function of f and f . (81) 

M 8 

A first-order correction to 77 opacity for a photon traveling at a small angle angle 9 S = arccos ^ s < 1 along the jet 
axis is obtained by replacing [i* with fj.*fj. s , which implicitly assumes that typical interact ions take place at a zimut h 
4>* = 7r/2. A detailed calculation of the 77 opacity from accretion disk photons is given bv lBecker fc Kafatod (|1995). 

3.2.3. Numerical Results 

The electron distribution is assumed to be well-represented by the iBand et al.l (|1993f )-tvpe function 

Nfo') = K' e H( 7 '; 7 ^ in ,^ ax ) {y-K exp(- 7 77 ') B"[(pa - PiH ~ 7'] + 

[(p 2 -PiK] P2 ~ Pl l'~ P2 eMP2~Pi)H[ 1 ' - (pa -pi) 7o ']}- (82) 

This distribution is essentially two smoothly joined power laws with power law number indices p\ and P2, and low- 
and high-energy cutoffs, j^in an d Imax? respectively, in the electron spectrum. For illustrative purposes, we take 
Pi = 2 and pi = 4, or a break by 2 units in the electron spectrum. This can be compared to a break by one unit 
expected for synchrotron and Thomson losses. Our approach is, however, to use the flaring synchrotron spectrum 
to imply the underlying flaring electron distribution without regard to specific acceleration and radiation processes. 
When the nonthermal electron distribution is obtained from analysis of blazar data, then the underlying jet physics 
that gives rise to the inferred electron spectrum can be examined. 

The total jet power in the stationary frame of the host galaxy is given by 

Pj = 2irR' 2 f3r 2 cu' tot = P j>par + P hB (83) 

(|Celotti fc Fabianlfl993l: iCelotti et al.ll2007t iFinke et ai1l2008fl . where u' tot is the total energy density in the jet, Pj lPar 
is the jet power from particles, and Pj.B is the jet power from the magnetic field. Here the factor of 2 takes into 
account that the jet is two-sided. For synchrotron-only emission, it is expected that the jet will be in equipartition 
and Pj, P ar ~ Pj,B, which will minimize the jet power. In order to explain the Compton-scattered component, however, 
the energy density in electrons can be different than the magnetic field energy density. 

We performed a parameter study by varying model parameters with respect to a baseline model, with baseline 
parameters given in Table 1. We consider a 10 8 Af Q supermassive black-hole jet source located at redshift z = 1. The 
jet is radiating at 10 3 i? g from the black hole and has Doppler factor 5^ = 25 and bulk Lorentz factor T = 25, so that 
the angle of the jet direction to the line of sight is 6 = 1/T. The mean magnetic field is 1 G, and the variability time 
t v = 10 4 s, corresponding to ~ lOx the light crossing time for the Schwarzschild radius of a 1O 8 M black hole. The 
jet opening angle 8j ss R' b /r < cfo£u/(l + z)r = 0.25 « 15° for standard parameters. While varying the parameters, 
the synchrotron spectrum was kept relatively constant by varying K' e and 7 ' following the relations given in ^3.2.11 
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Fig. 1. — Spectral energy distribution (SED) of accretion-disk/relativistic jet model using standard parameters for disk-jet system given 
by Table 1. Magnetic field is varied with the synchrotron flux remaining essentially constant, though with the peak synchrotron frequency 
varying between 2 X 10 and 10 eV. Separate spectral components (thin short dashed curves) are, from left to right, the synchrotron, 
accretion-disk, SSC, Compton-scattered accretion disk radiation, and second-order SSC components, for the high-magnetic field case, B = 4 
G. 

(except when cha nging angle). This was done with a \ 2 fitting technique to the baseline synchrotron spectrum (see 
iFinke et al.ll2008f ). The transition between the near field and far field takes place at f w T 4 , so that the baseline height 
of the jet is in the near field. 

The Compton-scattered accretion disk spectra are calculated from eq. (|70|) . Fig. [T] shows the effects of changing 
the magnetic field. With increasing B, fewer electrons are required to make the same synchrotron flux, so that both 
the first and second-order SSC flux, and the flux of the Compton-scattered accretion-disk component decrease with 
increasing B. In all of these models, the second-order SSC emission is overwelmed by the external Compton component 
and is not visible. The overall levels of all Compton-scattered components are oc B~ 2 . 

With increasing disk power fedd, the accretion-disk radiation and also the Compton-scattered accretion-disk compo- 
nent becomes progressively larger, as shown in Fig. [51 Note that the temperature of the disk photons at a given radius 

increases oc £^ d - Fig. [3] displays the dependence of the blazar SED on height f of the jet. As the blob's distance from 
the disk increases, the level of the Comptonized disk radiation decreases. In both the case of changing the fedd and 
r the SSC components are unaffected; eventually, the Compton-scattered disk flux falls below the SSC fluxes. The 
second-order SSC is visible when the external Compton flux decreases enough. Fig. [4] shows how variations in j m i n 
affect the low energy part of the synchrotron, SSC, and external Compton components. 

In Fig. [51 the Doppler factor is varied from the baseline value of 25, with the same approximate synchrotron flux. As 
<5d increases, the SSC component decreases, while the Compton-scattered accretion disk component increases. This 
behavior follows from eqs. (fl4|) and (f75|) . for the following reasons: For larger Doppler factors and fixed variability 
times, the radius becomes larger and fewer electrons are required to make the synchrotron flux at the same flux level, 
so the electron and internal photon energy densities decrease. Consequently the SSC flux decreases, but the larger 
value of Sd means that the external target photon field becomes more intense. Thus the Compton-scattered accretion 
disk radiation increases with increasing (5d- Also, as the Doppler factor increases, the SSC component is shifted to 
lower energies while the Compton-disk component is shifted to higher energies, in accordance with eqs. (|15[) and (|7S)) . 
respectively. 

Fig. [6] illustrates how the blazar SED is affected by changes in viewing angle. Contrary to Figs. - (O, we let the 
synchrotron flux vary. In this calculation, T — 25, and the observation angle and therefore Sd changes. It is interesting 
to note that the SSC flux varies least with changes in viewing angle, while the Comptonized disk component changes 
most rapidly. The ratios of the various components are explained by noting that Sd/T = 1,2/5,1/5, and 2/17 for 
6T = 1,2,3, and 4, respectively. The beaming factor for synchrotron radiation is oc S^ a oc 5^ in the flat portion of 
the vF u SED, with the convention that flux density F v oc v~ a . The relative magnitudes of the SSC components are 
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Fig. 3. — Same as Fig.fl] except that the jet height is varied. Separate spectral components are shown for the r = 3.2 X 10 4 R„ case. 

explained from eq. (fT"4|) . noting that / s it" oc Sfy, so that /'ssc oc 8^. The peak flux of the scattered disk component 
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Fig. 4. — Same as Fig. [T] except that the minimum electron Lorentz factor 7^, in is changed. Separate spectral components are shown 
for the "fmin = 600 case. 




Fig. 5. — Same as Fig. [T] except that the Doppler factor is varied with the synchrotron component essentially held constant. Separate 
spectral components are shown for the <5d = 15 case. 

the near-field regime, f^ CNF oc <5q, from eq. (|75|) . 
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Fig. 6. — Same i 
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Fig. 7. — Opacity of a photon emitted at various jet heights f in units of R g , subject to 77 pair product ion attenuation with photons 
from a Shakura-Sunyaev accretion disk with parameters given in Table 1. The parameter £ is defined in eq. 1811 . 

A larger viewing angle or smaller Doppler factor implies a smaller size scale of the emitting region for a given 
variability time t v . The characteristic Thomson scattering depth is n' e aTR' b cx R' b ~ 2 for a constant comoving electron 
spectrum N' e {^f') (described in our calculations by eq. [82]) in a region of size R' b . Consequently, the suppression of 
the SSC component due to Doppler boosting is offset by the increased scattering depth for larger observing angles. 

Fig. [7] shows a calculation of the accretion-disk opacity r 77 (ei,f) for photons traveling along the jet axis, using eq. 
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TABLE 2 

Jet Powers for Models in Figs. Q]-[6] a 



B [G] 


$D 


'Jmin 


r r-i r\A'\ — 1 1 

L B [10^ erg s i ] 


t r-i n4 c i 1 1 

Lpar [W^ erg s ] 


t r-i nd-^i — In 

Ltot [10 4J erg s ] 


1 


25 


100 


6.6 


11 


11 


1 


15 


100 


.85 


36 


36 


1 


20 


100 


2.7 


19 


19 


1 


30 


100 


14 


7.5 


7.6 


1 


35 


100 


25 


5.3 


5.5 


0.5 


25 


100 


1.6 


27 


27 


2 


25 


100 


26 


4.8 


5.0 


1 


25 


100 


110 


2.0 


3.1 


1 


25 


30 


6.6 


18 


18 


1 


25 


300 


6.6 


5.7 


5.7 


1 


25 


600 


6.6 


3.0 


3.0 



a First model is baseline model. 

(|80|) and the parameters given in Table 1. Here we assume that the Shakura-Sunyaev disk reaches to the innermost 
stable orbit (Ri — 6) of a Schwarzschild black hole. At larger jet heights, the accretion-disk opacity declines and only 
increasingly higher energy photons are subject to attenuation due to the threshold condition. The effect of increasing 
£ is to increase the mean accretion-disk photon energy at a given radius, from eq. (|66p , and therefore lower the energy 
at which 7 rays can be attenuated. At r = 10 3 , only > 10 TeV photons are attenuated by the disk radiation field (with 
^Edd = 1). The t 77 corrections have not been included in the SED calculations, but are only important at v > 10 25 
Hz. 

In all of our models the jet is particle-dominated (Pj, par 3> Pj,b) and the magnetic field is below its equipartition 
value. The total jet power, also divided into magnetic field and particle powers, is shown in Table [2l Thi s is typical 
of res ults from modeling the observed X-ray and 7-ray emission of TeV blazars with the SSC component (|Finke et al.l 
l2008h . 

3.3. External Isotropic Radiation Field 

This is the case treated bv lGeorganopoulos et alJ (|2001f ). where jet electrons scatter a surrounding external isotropic 
radiation field u*(e*,fi*) = u*(e*)/47r. By integrating eq. l|2"5)) over the angle variables, recognizing that dfj,*d(j)* = 
2nd cos %/} for this geometry, one obtains 

^ iso =iS^ r^^ 1 r^ftfcrj, (84) 

where Fc{q,T e ) is given by Jones's formula (|Jonesl fl968). eq. ([7]). Because the scattering is taking place in the 



stationary frame 

with T e = 4e*7 and, as before, e s = (1 + z)e. The limits j m in and "f ma x are given by 



- /or\ 

r e (l - e s / 7 



and 



, mm = \e s (l + yj 1 + ±- ) (86) 
Imax = e *^ S H{f* - e a ) + l2H(e s - e*) . (87) 



For the case of the cosmic microwave background radiation (CMBR), 

= Ubb {e*, II.) = ^ "J , (88) 

47r exp(e*/B) — 1 

where = k-QT/m e c 2 is the dimensionless temperature of the blackbody radiation field, and T — 2.72(1 + z) K. 
Substituting eq. flSSJ) in eqs. and (02]) gives 

/ c,CMB B= 3me^g% - r a ^7 7- 2 ^(7/fo)F c(g ,r e ). (89) 



2d| A 3 7 * exp( e */6) - 1 j lmm 

Eq. (|8"9")l gives an accurate spectru m of radiation made when jet electrons with number spectrum Nl(y ') Compton- 
scatter blackbody photons (see, e.g. iTavecchio et aLll200Ct iDermer fc Atovarj|2002l : iBottcher et al.ll2008l ). 
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Fig. 8. — Idealized geometry of the scattering region, approximated as a spherically-symmetric shell of gas with radial Thomson depth 
rx between inner and outer radii Ri and R a , respectively, and a density gradient defined by index £, where n e (R) oc R^ . 

3.4. Scattered BLR Radiation Field 

The BLR is thought to consist of dense clouds with a specified covering factor that can be determined from 
AGN studi es. These clouds in tercept central-source radiation to produce the broad emission lines in broad- 
line AGN (|Kaspi fc Netzerl I1999D . The diffuse gas will also Thomson scatter the central source radiation. The 
scattered radiation prov ides an important source of target photons that jet electrons sca tter to 7-ray energies 
(iSikora et al.lll994lll997f). This radiation field also attenuates 7 rays produced within the BLR ( Blandfo rd fe Leyinso: 
19951 iLevinson fe B landfordlll995l; IBottcher fe Dermeij[l995l : iDonea fe Protherodl200l iLiu fe Baill2006b lReimeijl200" 
Sitarek fc Bednarekll2008h . 

Here we calculate the angular distribution of the Thomson-scattered radiation from a shell of gas with density 



n e (R) = n 



H(R; Ri, R ) 



(90) 



extending from inner radius Ri to outer radius R (see Fig. [5] the calculation of fluorescence atomic-line radiation 
differs by considering dense clouds with a volume filling factor) . The shell is assumed to be spherically symmetric with 

a power-law radial density distribution and radial Thomson depth tt = ctt f R ° dR n e (R). The Thomson-scattered 
spectral photon density can be estimated by noting that a fraction ss rn e (r)<jT of the central source radiation with 
ambient photon density n p ^(e„; r) = N p h{t*)l 47rr 2 c is scattered, giving a target scattered radiation field 



n e (r)a T N p h(e*) 



Anrc 



(91) 



when Ri < r < Rq. Here N p h(e*) = i(e*)/('TieC 2 e I |,) is the central source photon-production rate, assumed to radiate 
isotropically, and L(e*) is its spectral luminosity. 
A more accurate calculation of the Thomson-scattered photon density is obtained by integrating the expression 



n ph {e*;R) 



dV 



n(e*; R) 



(92) 



47rx 2 c 

over volume, where dV = R 2 d(f>dfidR and x 2 = R 2 + r 2 — 2rR cos 9 (Fig. [5]). Assuming that the photons are isotropically 
Thomson-scattered by an electron without change in energy, so that n(e»;i?) = N p h(e^)aTn e (R) / (^■nR 2 c), then 



Tlph (6* , 



L(e*)a T 
87rm e c 3 e* 



(/// 



dR 



n e {R) 



8irm e c 3 e*r 



R 



R 



R-r 



(93) 



from a spherically-symmetric electron density distribution (cf. iGouldl 119791 : IBottcher fc Dermerl 119951 for a time- 
dependent treatment). In the case of an isotropic, uniform surrounding medium, n e (R) — no, C — 0j Ri = and 
R — > 00 in eq. ([90|) . and eq. (|93|) gives n p h(e r ;r) — 3.324A ? p / l (e»)crTno/(87r7-c), noting that the integral J °° du In |(1 + 
m)/(1 — u)\/u = 3.324. Thus the approximation given by eq. (f9Tj) is accurate to within a factor of w 2 for this case. 

The angle-dependent scattered photon distrib ution n n h(e*i M*i f) can be derived by imposing a (5-function constraint 
on the angle 6*» = arccos /x* in eq. ([92f (see also IDonea fc Protheroe]|2003f ) , so that 



n ph (e^.,^;r) 



%-Kcr 



du I dg 



g 2 + l- 



(94) 
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Fig. 9. — Angle-dependent density of Thomson-scattered radiation at height r along the jet axis. The radial Thomson depth of the 
spherically-symmetric scattering medium is tt = 0.01, and the scattering shell is assumed to extend from 10 2 to 10 5 Schwarzschild radii 
for a 10 8 Mq black hole. The solid and dashed curves show results with £ = and £ = —2, respectively. The values of r for are labeled on 
the dashed curves. The ratio of the luminosity of the photon source to the dimensionl ess p hoton energy is L»/e* = 10 44 ergs s —1 , and the 
calculation assumes that e* <g 1, so that scattering is in the Thomson regime (see eq. 1990 . 

after changing variables to g — R/r. From Fig. [5J m ax = 7r — f?», so that \i m in = cos 6 m ax — n — 0* and — //* < fx < 1. 
The law of sines gives R 2 (l — [i 2 ) = x 2 (l — with the result 

M^)= y f (1 ~^ . (95) 
y/1 + g 2 - 2 5A i 

Transforming the <5-function in fx* to a (5-function in g gives, after solving, 



where 



(units of AT are 1/L 3 ), and 



n p h(e*,^;r) = aT ^ ph e * jV(/x»,r) , (96) 
Hirer 

M{H*,r) = N[^,n{r)] = f dp n e (gr) ^1±°—^ (97) 

J -a, 5(1 — M ) 



g = = 2 2 ■ 98 

Fig. [5] shows the angle dependence of the scattered radiation field in the stationary frame for this idealized geometry 
for the parameters given in the figure caption. As can be seen, the scattered radiation field is nearly isotropic when 
r < Ri, and starts to display increasing asymmetry peaked in the outward direction at increasingly greater heights. 
When r > R , all scattered radiation is outwardly directed. Most of the scattering material is near the inner edge 
when C = —2, so that the intensity of the scattered radiation field is largest at r < E4. By contrast, the intensity of 
the scattered radiation field is not so enhanced towards the inner regions when £ = 0. 
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The Compton-scattered radiation spectrum is given, in general, by eq. (|34[) . Substituting eq. (]96|) for the angular 
distribution of the target photon source gives, for a monochromatic photon source 

N ph (e*) = L ° Sie * ~ e * o) , (99) 
m e cre* 



the vF v flux 

f EC, scat f \ _ ( 

127T 2 d 2 L r \e*J i max (-i,i_2e s / e ,) J^ ow V 



In this expression, 



7ic 



2 




(101) 



e*e s (l - 

(compare eq. |29j). Substituting eq. (|96p into eq. (f20|) for a monochromatic photon source, eq. (|99|) . gives 

r 77 (ei,r) = 3— / — / dfi* (1 - /x*) 7V(^», x)cr 77 (s) (102) 

o7rm e c e* j r a; j_i 

for the opacity of a photon with measured energy ei emitted outward along the jet axis at height r. The 77 opacity 
vanishes when e\ < l/[e*(l + z)\ due to the 77 pair-production threshold. 

3.4.1. Thorns on- Scattered Isotropic Monochromatic Radiation Field 

Substituting u*(e*,Q«) = u*o<5(e* — e*o)/47r for an isotropic monochromatic radiation field in eq. (f35|) gives, with 
5 — > 2 for the Thomson regime away from the endpoints of the spectrum, 

f T,*so = c*r e 6 u . / ^ / ^ t -2^ (7/( j d) (103) 

47Tfl L Jmax(-l,l-2e s /e») J^e s /2e t ) 



We consider Compton up-scattering, i.e., e s > e* for the power-law electron spectrum given by eq. (|47|) . The Thom- 
son approximation is only valid far from the endpoints. The asymptotes for the Thomson-scattered spectrum of a 
surrounding isotropic, monochromatic radiation field therefore becomes 



2c7rr 2 tt*o 



4ndl(p+ 1) " V e * 



and 



2 



9P+3 C7rr 2i3+P / n(3- P )/2 

/ e T ' 4S ° = : Ar, I if ' ( — ) , for 4e s (^7i) 2 < ^ « ^iltf ■ (105) 



(p+l)(p + 3) 47rrf 



Note the diffe r ent b eaming factors in the two asymptotes. Eq. (|105[) agrees with the Thomson expression derived by 
iDermer et al.l ()1997[ ). eq. (22), to within factors of order unity. 

3.4.2. Numerical Calculations 

We now present calculations of the SED of FSRQ blazars, including the external Compton-scattering component 
formed by jet electrons that scatter target photons which themselves were previously scattered by BLR material. This 
EC BLR component is not found in conventional synchrotron/SSC models of blazars, and so distinguishes standard 
model blazar BL Lacs and FSRQs. The inclusion of the 77 opacity through the scattered radiation makes the 
calculation self-consistent. 

For the simplified shell geometry depicted in Fig. [5J we calculate the external Compton scattering component using 
eq. (|100[) . and calculate the 77 opacity using eq. (]102|) . To simplify the calculations, the spectrum of the radiation 
scattered by the BLR is assumed to be monochromatic with energy w 50 eV, corresponding to the mean energy 
from the accretion-disk radiation for the standard model. Results of such a calculation are shown in Fig. [121 using 
parameters for a standard FSRQ blazar model given in Table 1. The constant-density BLR is confined between 100 
and 200 R g , and the Thomson depth through the BLR is tt = 0.1. In Fig. [lO^a), the emission region of the jet is 
far outside the BLR, at 1000 R g . Therefore most of the BLR photons encounter the jet from behind, and the rate 
of tail-on scatterings is suppressed by the rate factor. Consequently, the flux of the EC BLR component is much less 
that the flux of the EC Disk component. Fig. flOFb) presents a similar calculation, except that the BLR now extends 
to 10 3 R g , the same radius where the radiating jet is assumed to be located. As expected, the EC BLR component is 
significantly enhanced compared to Fig. [TOT a). 

There is very little 77 absorption by the accretion-disk radiation or the BLR radiation when the jet is found outside 
the BLR. On the other hand, when the jet is within the BLR, there can be significant 77 opacity. This is shown in 
Fig- [HI where we use the standard parameters for the jet and a BLR with tt = 0.01, except that the outer radius, R a , 
of the BLR is varied. The effects of 77 attenuation by the scattered radiation field are shown. When R a << r, the 
jet height, then the threshold is suppressed except for the highest-energy 7 rays. When the blob lies within the BLR, 
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Standard model FSRQ blazar with standard parameters given in Table 1. The optical depth through the BLR is tt = 0.1, 
ent in the density of the scattering material. The radiating jet is located at 10 3 ii 9 . Separate synchrotron, accretion-disk, 
SSC, EC Disk, and EC BLR components are shown, (a) The BLR extends from 100 to 200R 9 . (b) The BLR extends from 100 to 1000i? 9 . 




Fig. 11. — Model blazar SED comprised of synchrotron, accretion-disk, SSC, EC disk, and EC BLR components, including 77 attenuation 
by the scattered BLR radiation. The jet height r = W 3 R g . The outer radius, R , of the BLR is varied, keeping its Thomson depth, 
tt = 0.01, constant. 

r < R , the opacity from the scattered BLR radiation is significant and the 77 opacity large for 7 rays with e ~ 2/e*, 
which for our monochromatic radiation field with e* = 10~ 4 , is at v ~ 5 x 10 24 Hz (w 40 GeV). When R ^> r, the 77 
opacity, proportional to rn sc (e t ;r), declines oc R^ 1 for constant Thomson depth, as can be seen from eq. ([§1]) . Note 
that the EC 7 rays are not very sensitive to the changes in the BLR parameters, since most of the emission is formed 
by the EC disk component. 

The effect of changing the radial gradient of BLR scattering material is shown in Fig. [121 For a steeper density 
gradient, £ = — 2, and a constant Thomson depth, the material is concentrated near the inner edge of the BLR at Ri, 
so the changes in the 77 opacity are most dramatic when R Q « r due to geometric effects. When Ri <C r <?C R , the 
radiation field is essentially unchanged for different values of R , and so also is the 77 opacity. 
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Fig. 12. — Same as Fig. 1111 except that the gradient £ = —2 in the density distribution. 



We have presented accurate expressions for modeling synchrotron and Compton-scattered radiation from the jets of 
AGN that include target radia tion fields from th e accretion disk and BLR. This extends our technique for modeling 
synchrotron and SSC emission (|Finke et al.ll2008l) to include external Compton scattering in a relativistic jet of thermal 
radiation from the accretion disk, and accretion-disk radiation Thomson-scattered by electrons in the BLR. These 
formulae use the full Compton cross section and are accurate throughout the Thomson and Klein- Nishina regimes at 
any angle with respect to the jet axis, so can also be used to model 7-ray emission from radio galaxies, e.g., M87 
(jAharonian et al.l [2006). We also derive expressions for the 77 opacity through the same scattered radiation field 
that serves as a target photon source for the jet electrons. The expressions, eqs. (|80| and (|102p . for opacity from the 
accretion-disk and scattered radiation field assume, however, that the photon travels along the jet axis, though it is 
straightforward to derive the more general case. 

In the results presented here (Figs. [1] - [6] and [10] ~ I12|) we have chosen parameters for demonstration purposes 
that give an exaggerated EC component, particularly a high ^Edd- Lowering the disk luminosity lowers the radiation 
considerably, as seen in Fig. [5] The disk radiatio n field in FSRQ b lazars can be seen when the nonthermal blazar 
radiation is i n a low state, as in the cas es of 3C 279 (|Pian et al.lll999f) . 3C 454.3 (|Raiteri et alj [2007) and, most clearly, 
3C 273 (e.g.. lvon Montigny et al.lll997h . These observations can be used to assign the accretion-disk luminosity when 
modeling a specific blazar, though the disk brightness could also vary during the flaring epoch. 

The models presented here do, however, have limitations. They do not yet include enhancements from secondary 
cascade radiation initiated by pairs formed by 7 rays interacting with lower energy radiation from the disk and the 
BLR. For the parameters considered here, this would not make a significant difference in the calculated SEDs because 
the energy flux of the attenuated radiation is a small fraction of the escaping flux. But even in this case, the reinjected 
pairs from the attenuated radiation will be isotropized if the reinjection occurs outside the relativistic flow, and will 
then make only a small contribution to the Doppler-boosted radiation. 

Spectral features res ult from 77 abso r ption by BLR radiation, as seen in Fig. 1121 By assuming hard primary 7-ray 
emission components. jAharonian etafl |2008h argue that hard intrinsic spectra from blazars such as 1ES 1101-232 
(jAharonian et al.ll2007bh could be formed through 77 attenuation. If the primary TeV radiation originates from an 
underlying jetted electron distribution, then a consistent model requires that a 7-ray spectrum formed by Compton- 
scattering processes arises from the same radiation field responsible for 77 absorption. As our calculations show, soft 
Compton-scattered TeV spectra are formed due to Klein-Nishina effects in scattering, whether from the accretion 
disk or from p hotons scattered by the BLR. Thus either an extremely bright GeV component would be expected in 
the scenario of lAharonian et all (|2008[ ) (cf. Figs. 1111 and [T2"|), which would easily be detected with Fermi, or a hard 



7-ray FSRQ Blazar Analysis 



21 



primary 7-ray spectrum must originate from other processes. Moreover, if this explanation was correct, then blazars 
with stronger broad emission lines should have harder VHE 7-ray spectra than weak-lined BL Lacs. 

Cascade radiation induced by ultra-relativistic hadrons could inject high-energy leptons to form a hard radiation 
component, tho ugh a sufficiently dense t arget field for efficient photohadronic losses will itself severely attenuate the 
TeV radiation (| Atovan fc Dermerl [2003) . Depending on the underlying acceleration model, a distinctive hadronic 
signature could appear at 7-ray energies only, though one would expect that both electrons and protons would be 
accelerated synchronously. 

Our calculations also do not as yet include absorption by the diffuse extragalactic background light (EBL), which 
would be important for EGRET 7-ray FSRQs, whic h have a broad redshif t distribution with a mean value (z) ~ 1, 
larger than the mean value (z) ~ 0.3 for BL Lacs (|Mukheriee et al.l [l997). In cases where the Compton-scattered 
spectra decline steeply due to Klcin-Nishina effects (e.g., Fig. [J), the issue of EBL abso rption is secondary. In the on e 
FSRQ, 3C 279 at z = 0.536, detected from 80 to > 300 GeV energies with MAGIC (|MAGIC Collaboration l2008h . 
EBL effects cannot be neglected. For the models studied here, the soft calculated 7-ray spectra cannot in any case 
account for the measured, let alone intrinsic, spectrum of 3C 279. Based on the simultaneous optical - X-ray - VHE 
7-ray SED of 3C279 during the MAGIC detection, Bottcher, Marscher & Reimer (2008, in preparation) show that 
one-zone leptonic models have severe problems explaining the flare, and require either extremely high Doppler factors 
or magnetic fields well below equipartition. Correlated X-ray, Fermi, and TeV campaigns will offer opportunities to 
apply the results developed in this paper. 

In conclusion, we have derived expressions to model FSRQ blazars that self-consistently include 77 attenuation on 
the same target photons that are Compton-scattered by the relativistic electrons. By assuming that the lower energy, 
radio - UV emission is nonthermal synchrotron radiation, then the underlying electron distribution can be determined 
given the Doppler factor, magnetic field, and size scale of the emission region. The equations derived here can be 
used to calculate the 7-ray emission spectrum from SSC and external Compton scattering processes from this electron 
distribution. Separately, one can determine whether the inferred electron distribution can be derived from specific 
acceleration scenarios. Complementary to the technique of injecting electron spectra and cooling, this method can be 
applied to multiwavelength data sets to analyze high-energy processes in the jets of AGN. 

We thank the referee for asking us to address the potential importance of higher-order SSC fluxes, and for other 
helpful comments. The work of J.D.F. is supported by NASA Swift Guest Investigator Grant DPR-NNG05ED411 and 
NASA GLAST Science Investigation DPR-S-1563-Y, which also supported summer research by H.K. at NRL, and a 
visit by M.B. to NRL. C.D.D. is supported by the Office of Naval Research. 
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